How many CNOT gates does it take to generate a three-qubit state ? 
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The number of two-qubit gates required to transform deterministically a three-qubit pure quantum 
state into another is discussed. We show that any state can be prepared from a product state 
using at most three CNOT gates, and that, starting from the GHZ state, only two suffice. As a 
consequence, any three-qubit state can be transformed into any other using at most four CNOT 
gates. Generalizations to other two-qubit gates are also discussed. 

PACS numbers: 03.67.-a, 03.67.Ac, 03.67.Bg 

Quantum information and computation (see e.g. [H) is 
usually described using qubits as elementary units of in- 
formation which are manipulated through quantum oper- 
ators. In most practical implementations, such operators 
have to be realized as sequences of local transformations 
acting on a few qubits at a time. Whereas one-qubit gates 
alone cannot create entanglement, it has been shown that 
together with two-qubit gates they can form universal 
sets, from which the set of all unitary transformations of 
any number of qubits can be generated 0] ■ The complex- 
ity of a quantum algorithm is usually measured by assess- 
ing the number of elementary gates needed to perform 
the computation. The Controlled-Not (CNOT) gate, a 
two-qubit gate whose action can be written |00) — > |00), 
|01) ->■ |01), 1 10) -> |11), |11) ->■ 1 10) , is one of the most 
widely used both for theory and implementations. It can 
be shown that the CNOT gate together with one-qubit 
gates is a universal set Experimental implementa- 
tions of a CNOT gate (or the equivalent controlled phase- 



flip) have been recently reported using e.g. atom-photon 
interaction in cavities [H , linear optics 0] , superconduct- 
ing qubits 0, [|| or ion traps 0, Q. While large size 
quantum computers are still far away small platforms of 
a few qubits exist or can be envisioned in the framework 
of these existing experimental techniques. In most such 
implementations, two-qubit gates such as the CNOT are 
much more demanding that one-qubit gates. 

Theoretical quantum computation has been usually fo- 
cused on assessing the number of elementary gates to 
build a given unitary operator performing a given com- 
putation. Some works have tried to focus on two-qubit 
gates and to minimize their number in order to build a 
given unitary transformation for several qubits fioj ] . 
Still, unitary transformations in many applications are 
a tool to transform an initial state to a given state. It 
seems therefore natural to try and assess how costly this 
process is in itself. In this paper, we thus study the min- 
imal number of two-qubit gates needed to change a given 
quantum state to obtain another one. Of course, this 
number is necessarily upper bounded by the number re- 
quired for a general unitary transformation. We focus 



on the case of two and three qubits. For two qubits, we 
generalize the result of [ll[ and show that one CNOT 
is enough to go from any given pure state to any other 
one. For three qubits, we show that three CNOTs are 
enough to go from |000) to any other pure state, and 
that two CNOTs suffice if one starts from the GHZ state 
(|000) + |111))/V2. A corollary of the latter is thus that 
four CNOTs arc enough to go from any pure state to any 
other pure state. The number of CNOT gates required to 
go from a state to another defines a discrete distance on 
the Hilbert space. Given any fixed state the Hilbert 
space can be partitioned according to the distance to \ip). 
It is known that if stochastic one-qubit operations are 
used, entanglement of three [l2j and four [13| qubits fall 
into respectively two and nine different classes. Our clas- 
sification according to the number of CNOTs is different, 
although there are some relations. Our results generalize 
to other universal two-qubit gates, in particular to the 
iSWAP gate which has been shown to be implementable 
for superconducting qubits 14 1 . 

We consider pure states belonging to the 2 n - 
dimensional Hilbert space C 2 . The space of normal- 
ized quantum states is the sphere S 2 _1 . As the 
cost of one-qubit gates is negligible, we are interested 
only in equivalence classes of states modulo local uni- 
tary transformations (LU). We thus consider the sets 
£„ = 5 2 " +1_1 /C/(2)™ of states nonequivalent under LU. 
In the case of two and three qubits the dimension of 
£i and £3 was determined i n KEl Hfjj and their topol- 
ogy has been described in [17|. Throughout the pa- 
per we will make use of the one-qubit LU operations 



f\0 = exp(-i^f 
matrices acting on qubit k. 
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ation Ry k \t;) corresponds to 



In particular, the oper- 
a rotation of the qubit 
cos(p)|0) +shu»|l) i-> cos(<^ + £)|0) + sin(p + £)|1). 

Two-qubit states. Let us first consider the two-qubit 
case. It has been shown in [ll| that one can transform 
an arbitrary two-qubit state to 1 00) by using only 
one CNOT. Here we prove that the same holds for two 
general two-qubit states \ip) and \tp'). 

Proof: Since £2 is homcomorphic to [0,1], only one 



2 



parameter (e.g., one Schmidt coefficient) characterizes 
a state up to LU. More precisely, by LU each two- 
qubit state can be brought to the canonical form = 
cos^?|00) + sin(^|ll), which is just Schmidt decomposi- 
tion. We want to transform state with parameter ip 
to state with parameter ip'. When ip' ^ ip, we need 
at least one CNOT. It turns out that one CNOT is in 
fact sufficient, as can be easily seen by checking that the 



relation 
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(i) 



(<p'M) = W) holds (by 



convention, qubit 1 is the leftmost one). □ 

In contrast, one needs three CNOTs in general to con- 
struct a specific two-qubit unitary transformation Q. 
Transforming one state to another is thus clearly easier. 

We now turn to the three-qubit case. 

Classification with respect to 1 000) . We start with 
the case where we want to prepare a state \ip) from 1 000) . 
The distance (in number of CNOTs) from |000) to \tp) is 
a criterion for the difficulty to prepare \tp). We will show 
that this distance partitions the Hilbcrt space into four 
classes, and that any state can be prepared from |000) 
using three CNOT gates or less. We will examine each 
of these four classes in turn. 

Class 0: One needs zero CNOT gates to transform 
to 1 000) iff the state is of the product form \tp) = \aftj), 
where |a), \ft), I7) are normalized single qubit states (this 
is trivial, since only LU are used). 

Class 1: One needs one CNOT iff the state is of the 
form \ip) = |a)i|x)23 (i.e., it is bi-separable), where 
is an arbitrary entangled state of the last two qubits \li 

Proof: By LU can be transformed into canonical 
formic) h = |0)(cos<£|00)+sinv?|ll)). Applying CNOT 23 
to this canonical form we obtain |0)(cos c,c|00) +sin93|10)), 
i.e., state |0)(cos(^|0) + sin^l))^) which is in class 0. 
We can therefore reach |000) in a single CNOT step. 
Conversely, applying one CNOT gate on a state from 
class the state reached is bi-separable and therefore all 
states that need 1 step to get from |000) are of the above 
form. □ 

In their canonical form states from class 1 can be 
parametrized by a single real parameter <p. 

Class 2: One needs two CNOT gates iff the state 
is of the form \ip) = cos ip\aft^} + sin ip\a± ft' 7') , with 
(a\a±) = 0, \(ft\ft')\ < 1 and |( 7 | 7 ')| < 1 (if \(ft\ft')\ or 
I (7 1 7') I are equal to 1 then l^) belongs to class 1). 

Proof: We first bring the state by LU to the canon- 
ical form \ip) = cosi^lOOO) + sinp\lftj) . If the phases 
are absorbed into the definition of local bases, I7) 
can be written as cos£|0) + sin£|l). The rotation 
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followed by a CNOT13 yields a state 

(cos<y9|00) +smip\lft})\^'}, which is a state of class 1 from 
which we can reach state |000) in a single step. 

To prove the converse, we have to show that by using 
two CNOT gates one can reach only states of the form 
= cosip\aft^f) + smip\a±ft'j'), or states in class or 
class 1. Starting from class 0, we are in class 1 after one 



step. Any class 1 state can be written as |a)(cos <p|07) + 
c i? singly)). Applying CN0T 23 or CNOT 32 we get a 
state in class or class 1. Applying CNOT 2 i on the 
other hand we get cosip\aQj) + c 1 ^ smip\al~f') , where 
\a) = o~ x \a), which is indeed of the canonical form of 
class 2 states. Last possibility is applying CNOTi 2 . In 
this case it is better to write our state in the basis |±) = 
(|0) ± \1))/V2 for the second qubit. That is, any class 1 
state can be written as |a)(cos <p\ — 7) + e 1 ^ sint^l + 7')). 
Writing |a) = cos<//|0) + siny/|l), we then get after 
applying CNOTi 2 state cos <p(cos <p'|0) — c 1 ^ sin<//|l))| — 



7) + sin (pe 1 ^ (cos <p'\0)- 



-e 1 ^ smtp' 



l))|+7'), which is again 



of the canonical form expected. For CNOT 13 or CNOT 31 
the argument is similar. □ 
One needs 3 real parameters to describe states of class 
2 in their canonical form. Note that class 2 states con- 
stitute a subset of GHZ type states which are of (unnor- 
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malized) form \aftj) + \a' ft'^') 

Class 3: One needs three CNOT gates iff a state is 
not in class 0, 1 or 2. 

Proof: States not in the previous classes are of two 
types: (i) W-like states (according to the classifica- 
tion in 12j) for which the range of the reduced den- 
sity matrix of qubits 2 and 3 contains only one prod- 
uct state. Such states are of the (unnormalized) form 
\tp) = \aft-y) + |a')|x)23, where \x)23 is entangled and or- 
thogonal to \ftj). Under LU they can be written in the 
following canonical form [l2j 



LU 



cos<^|000) + sin^|a)(cosv?'|10) + sm<//|01)), (1) 



with |a) = cos£|0) + sin£|l) and, (ii) GHZ-like states 
with the canonical form 



|V) L - a 1 000) 



(2) 



where \a), \ft) and I7) arc real single qubit states 
parametrized by one parameter each and a, b are real 
parameters, one of which is fixed by normalization [l9j |. 
To exclude class 2 states we must demand that none of 
\a), \ft) and I7) be equal to |1). To exclude class 1 and 
states in ([I]),© p 23 must be of rank 2. W-like states ([1]) 
require 3 parameters. GHZ-like (J2J) states need 5, and 
thus are the generic states. 

First we show that by using single CNOT one can 
transform class 3 states to class 2. For W-like states ([T]) 
we just have to apply CNOT 23 to the canonical form (p} 
and we immediately get cost^OOO) + sin </?|a)(cos </?'|l) + 
sin<^'|0))|l) which is of class 2 (states on the third qubit 
arc orthogonal). For GHZ-like states ([2j) it is a bit more 
work. Note that GHZ-type states can be, by rearranging 
terms and after LU, written as cos ip\ 000) + siny)|l)|x)23 
(expanding |0) on the first qubit in ([2]) into |a) and \a±) 
and adding \a) part to the second term), where |x) 23 can 
in turn be expanded as |x) 23 = cos>p'\05) + sin <//|l<5')) 
with |(<^|5')l < 1 (otherwise state would be in class 
2). Finally, rotating third qubit brings the state to 
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cos<p|007'0+ sm< p|l)( cos ¥>'|00)+e 1? sin y>' 1 17') ) with real 
I7') = cos</|0) + sin g' | 1). After application of the ro- 
tation Rf ] f f - ^ followed by CNOT 23 we arrive at 

cos (^1007") + sin^|l)(cosv'|0) + smip'\l))\f) which 
is of class 2. 

Since the canonical forms of classes 0, 1, 2 and 3 span 
the whole Hilbert space and since the forms of classes 0, 
1 and 2 are the only ones that can be reached in 2 steps, 
it immediately follows that the states of canonical forms 
(H|) and require exactly three steps. □ 

Thus any state is at a distance less than or equal to 3 
from 1 000). 

Classification with respect to GHZ state. Let 

us now examine the distance to the GHZ state |GHZ) = 
(|000) + 1 1 1 1} ) /y/2. It turns out that any state is at a 
distance less than or equal to 2 from GHZ. We use the fact 
that any state in £3 can be characterized by a set of six 
polynomial invariants [2l|, • Following [24] , we denote 
the first three invariants by Ii = trp 2 , 1 < i < 3, where 
Pi is the reduced density matrix of the ith qubit. Any 
three-qubit state is equivalent under LU to its canonical 
form [22j 

Ao|000) + Aie i0 |lOO) + A 2 |101) + A 3 |110) + A 4 |lll), (3) 

where the six parameters Ao, . . . , A4, p label the state in 
£3 and ^2 A| = 1. If we take < p < ir the parameters 
in the canonical form ([3]) are unique. 

Class 0: States at distance from GHZ are states 
whose canonical form ([3]) is GHZ (this is trivial). 

Class 1: States at distance 1 from GHZ are states 
whose canonical form ([3]) is (up to relabeling) 



to be the same as for GHZ. Since GHZ is symmetric un- 
der permutation of the qubits, one can consider that the 
CNOT gate applied is CNOT 23 (up to relabeling of the 
qubits before applying the CNOT). In this case the in- 
variant Ii = i is conserved. The reduced density matrix 
of the first qubit of \tp) therefore verifies trpf = 5 which 
in turn (since it is a 2 x 2 density matrix) implies that 
pi = gl. All states with this property can be written 
as \1>) = ^\0)\ X ) + TfWIXx), with ( X \X±) = 0. The 
canonical form of such states is precisely (jH). □ 
Class 2: All other states are at distance 2 from GHZ. 

Proof: We show that all other states are at distance 

1 from states of canonical form States of canonical 
form (HJ) are characterized by Ii = h . Therefore we have 
to prove that any state \ip) not in class or 1 can be trans- 
formed to a state verifying I± = | with only one CNOT 
gate. The reduced density matrix of the first qubit is a 

2 x 2-matrix p\ = ^ ^* \ ^ m ' s ^ invariant is 

thus given by h = A 2 + (1 - A) 2 + 2\B\ 2 . It is equal 
to I iff A = 5 and B = 0, which (since B is complex) 
yields three equations. One-parameter rotations of the 
qubits before applying CNOT yield free parameters. It 
turns out that a solution to the three equations always 
exists. Indeed, the state \ip) can be reduced either to 
the canonical form {T]) or to the canonical form ([2]). If 
|V') is in canonical form {T]), one can check that the state 
CNOTi 2 i4 1) (0i)j4 2) (0 a )lV'}. with 



1 



— |000)+A 2 |101) 



A 3 |110) + A 4 |lll), 



(4) 



tan 20i 
tan 20 2 



(5) 



cot 2 (p + cos 2£ 
sin 2£ 

sin (£ + 20i ) cos p' sin 2p 
sin (2£ + 20i) cos 2<p' sin 2 ip — cos 2 ip sin 20i 



(with A4 7^ I/a/2) i.e., states with I\ = 5. 

Proof: It is straightforward to check that the state 
(|4]) has the same invariants (i.e. is the same up 
to LU) as CNOT 23 i?y 2) (02)-Ry 3) (03)|GHZ) with 2 = 
i arcsin(A 2 v / 2) and 3 = | arccos(A 3 v / 2), and thus is at 
distance 1 from GHZ. Conversely, if a state \ip) is at dis- 
tance 1 from GHZ, then one of its invariants I\, 7 2 , I3 has 



is such that I\ = ^. Let us now suppose that \tp) is in 
canonical form @. The one-qubit states \a), \/3), I7) can 
be written in the form cos<y9|0)+sin</?|l) with parameters 
respectively 951, pi and <p 3 . Normalization imposes that 
a 2 + b 2 + 2a6cos£cos tpi cos pi cos P3 = 1. One can check 

that the state CNOT 12 i?i 1) (x)i? y 1) (0i)i? y 2) (0 2 )|V'), with 
parameters of the rotations given by 



tan 20i 



tan 2% 



(b 2 — b 4 ) cos 2pi — a 2 (a 2 — 1 — 26 2 cos 2 p\ cos2p 2 + 26 4 cos 2 p>i sin 2 (2(^i) sin 2 ip s ) 
b 2 sin(2</?i) (l — a 2 — b 2 + 2a 2 cos 2 (£2(1 — a 2 sin 2 p% + b 2 cos(2c/?i) sin 2 </? 3 )) 



^ b 2 sin (2<pi + 20i ) sin 2p 2 + 2a,b cos £ sin (<pi + 26i ) sin p 2 cos P3 

a 2 sin 20i + b 2 sin (2ipi + 20i) cos2<y9 2 + 2a6cos£ sin (pi + 20i) cos pi cos v? 3 



2a6sin£ sint^i smpi cos<p 3 (a 2 sin20i — b 2 sin [2pi + 20i)) 
2a sin pi sin pi (2ab 2 cos pi cos pi + b(a 2 + b 2 ) cos £ cos y 3 ) ' 

I 



verifies Ii = \. □ The results above show that the number of CNOTs 



4 



needed to transform a state to another state is much less 
than the number to produce all unitary transformations. 
Indeed, according to [l(|, one needs at least 14 CNOT 
gates to produce any three-qubit unitary transformation. 
We also note that the best available algorithm act ually 
does not saturate the bound, needing 20 CNOTs [10| . 
Our procedure is explicit, and selects a specific unitary 
transformation leading from one state to the other using 
less CNOTs than a general unitary transformation. 

Let us now discuss the applicability of these results to 
physical systems. In an experimental context, our re- 
sults can be used to construct any desired state from the 
initial state which is easiest to produce with a given sys- 
tem. In some experimental set-ups, two-qubit gates are 
built from a nearest neighbor interaction (for example 
in In this case, the common procedure is to use 

additional SWAP gates to transfer the states of two dis- 
tant qubits to nearest neighbors before performing the 
CNOTs. However, in our case, due to the symmetry of 
GHZ state, one can go from any state to the GHZ state 
in two CNOTs without the need of any quantum SWAP. 
Thus even if only nearest-neighbors CNOTs are available 
for three qubits on a line, still only 4 CNOTs are enough 
to go from any state to any other state. If one starts 
from 1 000) and one allows for relabeling of qubits in the 
final state, 3 CNOTs are still enough to go to any state, 
except for the GHZ-type class 3 states, Eq. @, which 
need an additional CNOT in this architecture. 

Any two-qubit gate can be expressed in terms of 
CNOTs and one-qubit gates. Thus our result will im- 
ply a bound in the number of two-qubit gates needed to 
go from one three-qubit state to another, for any other 
choice of universal two-qubit gate. We note that another 
popular two-qubit gate is the iSWAP which is natural 
for implementations corresponding to a XY-intcraction. 
As the iSWAP can be expressed in terms of one CNOT 
and one SWAP gate plus one-qubit gates [3], our re- 
sults apply directly to this particular gate provided the 
SWAP can be made classically: the number of iSWAPs 
needed to transform three qubits is then the same as for 
CNOTs. This in particular arises when the physical im- 
plementation allows for a coupling between any pair of 
qubits, as swapping two qubits is equivalent to relabeling 
the qubits for all subsequent gates by interchanging the 
role of the qubits. An important example is the case of 
superconducting qubits coupled to each other via cavity 
bus [l4l ]. one of the most promising recent developments, 
where the resonance can be tuned to couple any pair. 

In conclusion, we have shown that one needs only 
three CNOTs plus additional one-qubit gates to trans- 
form 1 000) to any pure three-qubit states. If one starts 
from the GHZ state, only two CNOTs are enough, and 
thus one needs only four CNOTs plus additional one- 
qubit gates to transform any initial pure three-qubit state 
to any other pure three-qubit states. An interesting open 
question is to find out whether four is the maximal dis- 



tance between two three-qubit states (it should be at 
least three since |000) is at distance 3 from class 3 states). 
It would be also interesting to know how these results 
translate to mixed states, or to pure states using SLOCC 
instead of LU. At last, generalizations to higher numbers 
of qubits may pave the way to better optimization of 
quantum algorithms, which are usually described as uni- 
tary operators but are sometimes just transformations of 
a given state into another one. 
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